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Abstract 

We consider QCD radiative corrections to vector-boson production in hadron 
collisions. We present the next-to-next-to-leading order (NNLO) result of the hard- 
collinear coefficient function for the all-order resummation of logarithmically-enhanced 
contributions at small transverse momenta. The coefficient function controls NNLO 
contributions in resummed calculations at full next-to-next-to-leading logarithmic 
accuracy. The same coefficient function is used in applications of the subtraction 
method to perform fully-exclusive perturbative calculations up to NNLO. 
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The transverse-momentum (g^) distribution of systems with high invariant mass M (Drell- 
Yan lepton pairs, vector boson(s), Higgs boson(s) and so forth) produced in hadronic collisions is 
important for physics studies within and beyond the Standard Model (SM). 

The computation of these distributions in perturbative QCD is complicated by the presence of 
large logarithmic contributions of the form ln(M 2 / q^) that need to be resummed to all perturbative 
orders in the QCD coupling as- The method to perform the resummation is known [U EJ El H] , 
including recent developments on the discovered and resummed effects |6] due to helicity and 
azimuthal correlations in gluon fusion subprocesses. The structure of the resummed calculation 
is organized in a process-independent form that is controlled by a set of perturbative functions 
with computable 'resummation coefficients'. All the resummation coefficients that are process 
independent are known since some time [7J El El EI] up to the second order in as, and the third- 
order coefficient has been obtained in Ref. [TT]. The complete computations of the second- 
order resummation coefficients have been carried out in Refs. [12] and [T3] for two benchmark 
processes, namely, the production of the SM Higgs boson through gluon fusion and vector boson 
production through the Drell-Yan (DY) mechanism of quark-antiquark annihilation. The explicit 
analytic expressions for the 0(a\) hard-collinear resummation coefficients in the case of SM Higgs 
boson production in the large-mfo P limit have been presented in Ref. [H]. This paper parallels 
Ref. [14]: we concentrate on single vector boson production, and we present the corresponding 
analytic expressions of the second-order hard-collinear coefficient functions H^. 

QCD predictions for vector boson production at hadron colliders are important for present and 
forthcoming studies at the Tevatron and the LHC. Resummed calculations of the qx spectrum of 
vector bosons and of related observables are presented in Refs. [U]-[2S]. Calculations for vector 
boson production at the fully-exclusive level with respect to the accompanying QCD radiation 
have been carried out in Refs. [271 (131 I2H] up to the next-to-next-to-leading order (NNLO) in 
perturbative QCD. 

In this paper we compute the hard-collinear coefficient function T^S 2 ' and, thus, the complete 
analytical expression of the NNLO cross section for vector boson production in the small-g^ 
region. These results have a twofold relevance, in the context of both resummed and fixed-order 
calculations. 

The knowledge of can be implemented in resummed calculations at full next-to-next-to- 
leading logarithmic (NNLL) order to achieve uniform NNLO accuracy in the small-gr region. In 
the case of vector boson production, this implementation has been carried out in Ref. [22J by using 
the impact-parameter space resummation formalism developed in Refs. [29, 30J. This formalism 
enforces a unitarity constraint and thus it guarantees that (upon inclusion of H^) the resummed 
qr spectrum returns the complete NNLO total cross section after integration over qT- 

The subtraction method of Ref. [12] exploits the knowledge of transverse-momentum resum- 
mation coefficients at 0(a^) to perform NNLO calculations at the fully-exclusive level. The 
Higgs boson coefficient functions presented in Ref. [H] were used in the numerical computations 
of Refs. [121 EI]- The coefficient functions presented in this paper are precisely those that are 
needed for the actual implementation of this subtraction method in DY-type processes: they are 
used in Refs. [13] and [32] for the NNLO numerical computations of vector boson production and 
of associated production of a Higgs boson and a W boson. The diphoton NNLO calculation of 
Ref. [33] also uses part of the results of the present paper to treat the quark-antiquark annihilation 
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subprocess qq — > 77. 



This paper is organized as follows. We first introduce our notation and illustrate the NNLO 
calculation of the vector boson cross section at small values of qx- Then we recall the transverse- 
momentum resummation formalism. Finally, we present our NNLO results in analytic form and 
the relation with the qx resummation coefficients at 0{a\). 

We briefly introduce the theoretical framework and our notation. We consider the production 
of a vector boson V (V = W ± ,Z and/or 7*) in hadron-hadron collisions. We use the narrow 
width approximation and we treat the vector boson as an on-shell particle with mass M. The 
QCD expression of the vector boson transverse-momentum cross sectionQ is 

da , , , . \ f 1 , f 1 , „ . . ,os „ , . da, 
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(g T ,M, S ) = V/ dzi / dz 2 f a/hl {z u M 2 )f b/h2 {z 2 ,M 2 ) ^{q T ,M,s = z 1 z 2 s;a s {M 2 )) 



(1) 

where f a /hi(x, Hf) (a = <?/, qf, g) are the parton densities of the colliding hadrons {h\ and h 2 ) at the 
factorization scale fj,p, and da ab /dq T are the partonic cross sections. The centre-of-mass energy 
of the two colliding hadrons is denoted by s, and s is the partonic centre-of-mass energy. We 
use parton densities as defined in the MS factorization scheme, and as(/^ij) is the QCD running 
coupling at the renormalization scale in the MS renormalization scheme. In Eq. (CQ) and 
throughout the paper, the arbitrary factorization and renormalization scales, fip and are set 
to be equal to the vector boson mass M. 

The partonic cross sections da ab jdq\ are computable in QCD perturbation theory as power 
series expansions in a${M 2 \ We are interested in the perturbative contributions that are large 
in the small-gr region (qx <C M) and, eventually, singular in the limit qx — >■ 0. To explicitly 
recall the perturbative structure of these enhanced terms at small qx, we follow Ref. [H] and we 
introduce the cumulative partonic cross sectiorJi] 

dq 2 ^(q T ,M,s = M 2 /z;a s (M 2 ))^ £ z a% R^ ab (z, M/Q ; a s (M 2 )) , (2) 

qT c=g f ,g f , 

where the overall normalization of the function R v is defined with respect to er^ t v , which is the 
Born level cross section for the quark-antiquark annihilation subprocess qf(jf — > V (the quark 
flavours / and /' are equal if V — Z, 7*). The partonic function R v has the following perturbative 
expansion 

00 

RZ^ ab (z, M/Q ; a s ) = 5 ca 5 5b 5(1 - z) + ^ f^)" R V J? ab {z, M/Q ) . (3) 

n=l 

The next-to-leading order (NLO) and NNLO contributions to the cumulative cross section in 
Eq. (T5]) are determined by the functions R v ^ and R v ^ 2 \ respectively. The small-gr region of the 



^If V — 7* or if the vector boson V is not an on-shell particle, the transverse- momentum cross section da/dq^ 
has to be replaced by the doubly-differential distribution M 2 da/dM 2 dq^ , where M is the invariant mass of V . 

•^In our notation, the subscripts c and c denote a quark and an antiquark (or viceversa) that do not necessarily 
have the same flavour. The flavour structure depends on the produced vector boson V and it is (implicitly) specified 
by the specific form of the Born level cross section u^y ■ 
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cross section da ab /dq T is probed by performing the limit Qq <C M in Eq. (j2J). In this limit, the 
NLO and NNLO functions R v ^ and R v ^ have the following behaviour: 

WQo) = « ££2*C*) + + J^C*) + 0(Ql/M 2 ) , (4) 

+ Ag^s) + ^(s) + 0(Q^/M 2 ) , (5) 

where l = lia(M 2 /Ql). In Eqs. (j4]) and (jSJ), the powers of the large logarithm l are produced 
by the singular (though, integrable) behaviour of da ab jdq\ at small values of q?. The coefficients 
^(i; m ) ( w ith m < 2) and i?( 2;m ) (with m < 4) of the large logarithms are independent of Qq] 
these coefficients depend on the partonic centre-of-mass energy s and, more precisely, they are 
functions of the energy fraction z = M 2 /s. As is well known (see also Eq. (J7|)), the logarithmic 
coefficients _R( n;m ) do not depend on the specific vector boson that is produced by qq annihilation 
and, therefore, we have removed the explicit superscript V (i.e., _R y ( n;m ) = _R(™; m )). 

In this paper we present the result of the computation of the cumulative cross section in 
Eq. (J2J) up to NNLO. The partonic calculation is performed in analytic form by neglecting terms 
of 0(Ql/ M 2 ) in the limit Qq <C M. Therefore, we determine the coefficient functions R^ n ' mS) (z) 
in Eqs. (jl) and ©. 

To perform our calculation, we follow the same method as used in Ref. [H] to evaluate the 
transverse-momentum cross section for Higgs boson production. The qx integration in Eq. (J2J) is 
thus rewritten in the following form: 

f Q2 °, 2 da ab f +ao 2 da ab A f + °° 2 da ab . 

/ «<?t TT s;a s ) = / dq T —^-{q T , M,s;a s ) - dq T — — (g T , M, s; a s ) 

Jo dq T J Q dq T J Q 2 dq T 

poo r+oo j*. 

= 5-( tot) (M,s;a s ) - / dq 2 T \ dy f (y, q T , M, s; a s ) , (6) 
JqI J-oo dy dq T 

where cr^ ot ^ is the vector boson total (i.e. integrated over q?) cross section and da ab /dy dq T is 
the corresponding doubly-differential cross section with respect to the transverse momentum and 
rapidity (y is the rapidity of V in the centre-of-mass frame of the two colliding partons a and b) 
of the vector boson. The total cross section a^ b (M, s; as) is known [51] in analytic form up to 
NNLO (i.e., up toO(a§4 0) )). In the region of large or, more precisely, non- vanishing values of q?, 
the differential distribution da ab / 'dy dq T is also known [35] [361 EI] i n analytic form up to 0(a\ay ). 
Using these known results and exploiting Eq. (jSJ), we can compute the cumulative partonic cross 
section up to the NNLO. Note that qT > Qq in the last term on the right-hand side of Eq. (j6J). 
Therefore the corresponding integration of the expression da ab / dy dq T [55] [56] [57] over y and q T 
is finite as long as Qq ^ 0: using the explicit expression of da ab / dy dq T frorr|§ Ref. [36], we carry 
out the integration in analytic from in the limit Qq <C M (i.e^ we neg lect terms of O(Q 2 /M 2 ) on 
the right-hand side of Eq. (jf)])). The result of our calculation.] confirms the logarithmic structure 

s We list some typos that we have found and corrected in some formulae of Ref. [3B]. In Eq. (2.12), Bf has 
to be replaced by J3f G + C| G , and Cf G has to be replaced by C| G . In Eq. (A. 4), two signs have to be changed: 
-B' G has to be replaced by —B qG , and A qG has to be replaced by —A qG . In the first line of Eq. (A. 10), the term 
C F (/« - Is- ft) has to be replaced by C A (f u -f 8 -f t ). 

^Some technical details related to the limit Q <§C M are illustrated in Ref. |14j . 
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in Eqs. @j and (J5J), and it allows us to determine the NLO and NNLO coefficients _R( 1;m ) (with 
m < 2) and _R( 2;m ) (with m < 4) of the cumulative cross section in Eq. (j2J). 

The results of the coefficient functions _R( n;m ) (z) are conveniently expressed in terms of transverse- 
momentum resummation coefficients. Therefore, before presenting the results, we recall how these 
functions are related to the perturbative coefficients of the transverse-momentum resummation 
formula for vector boson production [3]. This relation also shows that from the knowledge of 
Eq. ([5]) we can fully determine the NNLO rapidity distribution of the vector boson in the small-gr 
region. 

To present the transverse-momentum resummation formula, we first decompose the partonic 
cross section da a b/dq T in Eq. (GQ) in the form da a b = da^ b mg ^ + da^ b g \ The singular component, 
da^ b , contains all the contributions that are enhanced at small qr- These contributions are pro- 
portional to S(q T ) or to large logarithms of the type l/q T hi m (M 2 /q T ). The remaining component, 
da* b , of the partonic cross section is regular order-by-order in «s as qx — > 0: the integration of 
d&ab /d<lT over the range < qx < Qo leads to a result that, at each fixed order in a s , vanishes in 
the limit Qq — > 0. Therefore, da^ b g ^ only contributes to the terms of 0(Ql/M 2 ) on the right-hand 
side of Eqs. (j3J and ©. 

Inserting the decomposition da a b = da^ b + da^ b g ^ in Eq. (JTJ) , we obtain a corresponding 
decomposition, da = da^ smg ^ + da^ reg \ of the hadronic cross section. The transverse- momentum 
resummation formula for the singular component of the qr cross section at fixed value of the 
rapidity y (the rapidity is defined in the centre-of-mass frame of the two colliding hadrons) of the 
vector boson reads [3j H] 

—^(y, qT ,M,s) = — J2 a £v db-J (bq T )S q (M,b) 

y qT S c=q f ,q fl J ° 

X E / — / — i HFC ^]c- C - aia2 fa l/hl (Xl/ZlX/b 2 )fa 2/h2 (x 2 /z 2 ,bl/b 2 ) , (7) 
ai,oa 1 Jx * 2 

where the kinematical variables x^ (i — 1, 2) are x\ = e +y M/ ^Js and x 2 = e~ y M/ yfs. The integra- 
tion variable b is the impact parameter, Jo(&9t) is the Oth-order Bessel function, and 60 — 2e~ 7B 
(je = 0.5772... is the Euler number) is a numerical coefficient. The symbol \H F CiC-2\ c _ a ^ 
shortly denotes the following function of the longitudinal- momentum fractions Z\ and z 2 : 

[ HDYC i C *Ua ia2 = H r(MM 2 )) C cai ( Zl ^s(b 2 /b 2 )) C 5a2 (z 2] a s (bl/b 2 )) , (8) 

where H^ Y (as) and C ca (z;as) (c = qf,q~f) are perturbative functions of as (see Eqs. (TT2T) (TT3|) ) . 

The quark form factor S q (M,b) in Eq. (JT)) is a process-independent quantity [3], H]- Its 
functional dependence on M and b is controlled by two perturbative functions, which are usually 
denoted as A q (as) and B q (a$) (see, e.g., Ref. [6] that uses the same notation as in Eq. (JTJ)). Their 
corresponding n-th order perturbative coefficients are A q n) and B { q n) . The coefficients A q 1] , B { q l \ 
A q 2 ^ [7] and B q 2 ^ [8] are known: their knowledge fully determines the perturbative expression of 
S q (M, b) up to O(al). 

The perturbative function H q DY (as) in Eq. (jSJ) is process dependent, since it is directly related 
to the production mechanism of the vector boson through quark-antiquark annihilation. How- 
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ever, H® Y is independent of the specific type of vector boson V (V = W ± ,Z, 7*), and we have 
introduced the generic superscript DY. 

The partonic functions C q . a and C qfa in Eq. flB]) are instead process independent, as a conse- 
quence of the universality features of QCD collinear radiation. Owing to their process indepen- 



dence, these partonic functions fulfil the following relations: 

C qfqfl (z;a s ) = Cq f q fl (z;a s ) = C qq (z;a s ) 5 fr + C qq ,(z;a s ) (1 - %/) , (9) 

C qfq - f/ (z;oi s ) = Cq f g f/ (z;a s ) = C qq (z;a s ) 5 ff > + C qq >(z;a s ) (1 - 5 ff ) , (10) 

C qf9 (z;a s ) = Cg f g (z; a s ) = C qg (z;a s ) , (11) 



which are a consequence of charge conjugation invariance and flavour symmetry of QCD. The 
dependence of the matrix C ca on the parton labels is thus fully specified by the five independent 
quark functions C qq , C qq i, C qq , C qq i and C qg on the right-hand side of Eqs. (T9|) (TTTj) . 

We recall that the function if^ y (o;s), the quark functions C qa (as) and the perturbative func- 
tion B q (as) of the quark form factor are not separately computable in an unambiguous way. Indeed, 
these three functions are related by a renormalization-group symmetry [4] that follows from the 
6-space factorization structure of Eq. ([7]). The unambiguous definition of these three functions 
thus requires the specification^ of a resummation scheme [I]. Note, however, that considering the 
perturbative expansion^ of Eq. ([7]) (i.e., the perturbative expansion of the singular component of 
the qr cross section), the resummation-scheme dependence exactly cancels order- by-order in as- 

The perturbative expansion of the quark functions C qa (as) and of the vector boson function 
H q Y (as) is defined as follows: 

00 

C qa (z;as) = S qa 6(l-z) + J2{^) n C q n J(z) , (a = g,q,q,q',q') , (12) 

n=l 

00 

Os) = i + E(v) n < y(n) - (13) 

n=l 

The first-order coefficient function Cq~ g \z) is independent of the resummation scheme; its expres- 
sion is [8] 

C${z) = \z{\-z) . (14) 
The first-order coefficients C^l(z), Cq^(z) and C^,(z) vanish, 

while the coefficients C q )J (z) and H q Y fulfill the following relation [8j [10] : 

+ l H ? Y(1) *a -*) = t((t- 4 ) 6(1 - z ^ + l ~ z ) ■ ( 16 ) 

^The reader who is not interested in issues related to the specification of a resummation scheme can simply 
assume that (as) = 1 throughout this paper. The choice H® Y (as) — 1 is customarily used in most of the 
literature on qj< resummation for vector boson production. 

*The resummation-scheme dependence also cancels by consistently expanding Eq. ([JJ in terms of classes of 
resummed (leading, next-to- leading and so forth) logarithmic contributions [29] . 
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The separate determination of Cq q \z) and Hq Y ^ requires the specification of a resummation 
scheme. For instance, considering the resummation scheme in which the coefficient Hq Y ^ van- 
ishes, the right-hand side of Eq. (fl6l) gives the value of Cj^j (z), and the corresponding value of 
the quark form factor coefficient B q is explicitly computed in Ref. [8]. The computation of 
the second-order coefficients Cqq , C q 2 q ], , C^l, Cq 2 g and Hq Y ^ is the aim of the calculation 
described in this paper. 

To the purpose of presenting the NNLO results for the cumulative cross section in Eq. (T5]), we 
also define the following hard-collinear coefficient function: 



Ufg^-abi^ «s) = H q Y (a s ) \ dz l / dz 2 5{z - z x z 2 ) C qa (zi, a s ) Cg b {z 2 ; a s ) 

Jo Jo 



(17) 



which is directly related to the coefficient function in Eq. (jSJ). The function T-L DY depends only 
on the energy fraction z, and it arises after integration of the resummation formula ([7j) over the 
rapidity of the vector boson. Note that HP Y is independent of the resummation scheme [3]. The 
perturbative expansion of the function T-L DY directly follows from Eqs. ffT21) - ffT51) . We have 



oo 

nfql ab (z; a s ) = 5 qa 5q b 5(1 - z) + £ (^) B «J3 



n=l 

where the first-order and second-order contributions are 



= S qa Sq b 5(l -z)H° Y ^ +5 qa cf b \z) + 6q b C^(z) , (19) 



Cw = *«■ *** - *) < y(2) + °S (*) + c $ (*) 

In Eq. (1201) and in the following, the symbol £g> denotes the convolution integral (i.e., we define 
(g g> = j" Jg 1 cfz 2 S(z - z x z 2 ) g(zt) h(z 2 )). 

In the limit Qq <C M, the perturbative expansion of the cumulative partonic cross section in 
Eq. (|2J) can directly be related to the resummation coefficients of Eq. (j7j). We refer the reader to 
Ref. [H] for a concise illustration of this relation and to Ref. |29j for more technical details. The 
NLO and NNLO functions R v ^ and R v ^ in Eqs. fl3| and ([5]) have the following expressions: 

M /Qo) = ll CV) + V%S£\z) + K£2(z) + 0{Ql/M*) , (21) 



+ - 4C3 SSS 3) (*)) + 0(<K/J»O , (22) 

where we have used the same notation as in Ref. [29J. The explicit expressions of the coefficient 
functions S^j"^ (z) in terms of the resummation coefficients are given in Eqs. (63), (64), (66)- 
(69) of Ref. (we have to set hr = /ip = Q = M, where /ir, ^r and Q are the auxiliary scales 
of Ref. [2S]) and are not reported here. The coefficients "H^L^ an d W^Sab are exactly those in 
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Eqs. dll]) and (I2"0~|) (they are also given in Eqs. (65) and (70) of Ref. J2U]) The first-order terms 
^q^-ab an d ^^-06 depend on the quark form factor S q (M, b). The second-order terms S^L^, 

depend on H^J^b anc ^ 011 the Q uar k form factor S q (M, b) up to O(ag). The numerical coefficient 
~ 1.202 ... is the Riemann (^-function) on the right-hand side of Eq. (1221) originates from 
the Bessel transformations (see, e.g., Eqs. (B.18) and (B.30) in Appendix B of Ref. [29|). 

We now document our results of the NNLO computation of the cumulative partonic cross 
section. Using Eqs. (J2U) and (J22I), the results for _R y(1) and R v{2) allow us to extract and 
1lDY(n) U p ^ Q Q(^ a 2y The explicit result of the NLO function R v ^(z) confirms the expressions of 

^gf-oli ( z )i ^qX-ab \ z ) aric ^ ^^-06 '> as predicted by the resummation coefficients at C(as). 
At NNLO, the present knowledge [TJ [5] of the resummation coefficients at 0(a^) predicts 

the expressions of the terms ^^}^™\z), with m = 1,2,3,4. Our result for the NNLO function 
R v ( 2 \z) confirms this prediction, and it allows us to extract the explicit expression of the second- 
order coefficient function 
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We obtain 

v dy(2), s r r j fKs 101\( 1 \ /59C 3 1535 2157T 2 tt 4 \ 

l + z 2 f Li 3 (l-z) , T . , , Li 2 (z)log(*) U.,,, n , 1 



, _ : , , ■ Li 3 (^) - ZV 7 2 " V 7 - ^l-L-O log(l •:)••- log'i: 

- ± log 2 (l - *) log(*) + ^ Ml - *) - y) + 77^ ( - \ (11 - 3^ 2 ) C 3 

- ^ (-^ 2 + 12* + 11) log 2 (.) - 1 (83, 2 - 36, + 29) log(z) + ^ 

JU 2 (z) 1 . A ,+ 100 1 . A 

+ (!"*) (^^^ + g lQ g(l - z ) 1 °Z( Z )) + + ~ A Z lQ g(l " *) J 

+ H + J_ (192C 3 + 1143 - 1527T 2 ) 5(1 - *) 

+ ^(f5) log(2)(31ogW + 10) + To8 { ~ 19z ~ 37) } 

^ 2 fl / irA 511 67tt 2 17tt 4 \ 

1 + , 2 W l - z) 5Li 3 (,) + 1 _ + 3Li 2 (,) Iog(z) 



1 - 2 V 2 2 2 v/ov 7 2 

+ \ \og{z) log 2 (l - *) + \ \oi(z) log(l - ^) - ^tt 2 log(l - z) + ^ 

+ (1 - z) (-Li 2 (z) - | log(l - z) log(z) + ^ - ^ + 1 (1 + z) log 3 (z) 
+ y__7 Q (- 2 ^ + 2^ + 3) log 2 ^) + i (l7z 2 - 13z + 4) log(z)) - | log(l - 2) } 
+ C F {1(1 - z) (2, 2 - * + 2) + I log(l - s) log(s) - ^ 



' [l-z) (l36z 2 - U3z + 172) - — (8z 2 + 3z + 3) log 2 (z) 



(23) 



216,2 /v ' 48 

+ ^ (32z 2 - 30z + 21) log(z) + ^(1 + z) log 3 (z)| , 

= ~ 2) (2 ^ 2 " z + 2) ( Li2W + log(1 ~ z) log(z) ~ if) 

+ 432^ (1 " Z) (l3fe2 " ^ + 172) + i (1 + hg3{z) 

- 1 (8z 2 + 3^ + 3) log 2 ^) + 1 (32z 2 - 302 + 21) log(z)} , (24) 



8 



«SSw = c F (c, - l ,c A ) (i±^ f *m=£> + Li3(z) + Li3 ^ _ ^i-^m 



1 + Z \ 2 + z 

_ Li 2 (*)Mz) _ ^ log3(z) _ 1 log3(1 + z) + 1 lQg(1 + z) lQg2(z) 
+ £ Ml + *) - ^ V (1 - z) ( ^# + J log(l - z) log(s) + ^ 



12 ov 4 / v V 2 2 ov / ow ■ g 

- 1(1 + z) (!A 2 (-z) + log(z) log(l + z)) + ^ (z - 3) + 1(11* + 3) log(z) 

+ ^(1 - *) (2* 2 - * + 2) (Li 2 (*) + log(l - z) \og{z) - ^ 
+ - *) (l36z 2 - 143,2 + 172) - ^ (8z 2 + 3z + 3) log 2 ^) 



(25) 



+ 1 (32z 2 - 30z + 21) hg(z) + ^(l + z) \og 3 (z) | , 

Mq^ll'iz) = WqqJ^qi{z) , (26) 

n%®(z) = cu| - ^(i - *) (ii^ 2 - 2 + 2) Li 2 (i - z) 

+ (2^ 2 - 2, + 1) (iMlzf) _ l Li2 (i _ z) log(1 _ ,) + _L log 3 (1 _ ,) 



+ (2, 2 + 2, + l) ( + ^Mil - U ^~ Z l l0g(z) - 1 l 0g 3 (1 + z) 

K ' \ 8 4 8 24 

+ ^ log 2 ^) log(l + z) + ^tt 2 log(l + z^J + \z{\ + z)Li 2 (-z) + zU 3 (z) 

1 3 149* 2 

- -zU 2 (l - z) \og{z) - zU 2 {z) \og{z) - - (2z 2 + 1) Cs - 

- - (44z 2 - 12* + 3) log 2 (z) + - (68^ 2 + 6tt 2 z - 30* + 21) log(z) + — + 

+ itL + + 1} iog3(2) - r iog(i - 2) hg2(z) z)zhg2(i - 2) 

1 1 35 

+ -z(l + z) log(l + z) log(z) + —(3 - 4z)zlog(l - z) - — 

+ C F { (2. 2 -2z+l) ((a - - ^ + lLi 2 (l - *) log(l - z) 

+ Li2(g) g l0g(g) - 1 log 3 (l - z) + 1 log(z) log 2 (l - z) + 1 log 2 (.) log(l - *; 

3^ 2 1 1 

" T" " 4 " 2z + ^ log3{z) + el ( ~ 8 ^ + 122 + x) log2(2) 

+ - (-8z 2 + 23* + 8) log(z) + ^vr 2 (l - z ) z + -± + -(i- z )z log 2 (l - *) 

- 1(1 - z)zlog(l - ^) log(z) - ^(3 - 4z)*log(l - z) - A | , (27) 



where C F = (iV c 2 - 1)/(2JV C ), CU = iV c (JV C is the number of colours in SU(N C ) QCD), n F is the 
number of quark flavours and Lifc(z) (A; = 2, 3) are the usual polylogarithm functions, 

U 2 (z) = - / - ln(l-t) , Lis(z)= / - ln(t) ln(l - *f) . (29) 
Jo * Jo t 



We comment on the vector boson results in Eqs. ([23]) - ([28]) and on the ensuing determination 
of the second-order coefficients C { q 2 q \ Cf\, Cf), Cfi Cf) and Hq Y{2) in Eqs. 0X2} and (H3J. 
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The parton matrix %^L^2 is completely specified by the six entries^ in Eqs. 



quark-quark functions H^L^, H 
the gluon-gluon function Hqqt-gg- 



DY(2) 
~qq<-qq' 



the 

^■qqX-qqi WaZ-al'i tlie quark-gluon function HfZ-Zl and 



qq^qq 



"99-*-99 



Using Eq. (I2U]) . in the gluon-gluon channel we have 



njDY(2), x 

rL qq<-gg\ 4 ') 



(30) 



We see that the second-order coefficient function Tif^g^z) 



is fully determined by the resum- 



mation coefficients at C(as). Using the value of C^p in Eq. (JHj), the expression on the right-hand 



side of Eq. (I30p is in complete agreement with the result in Eq. (j28j) . Therefore, our explicit com- 
putation of the NNLO partonic function R^^} gg represents a consistency check of the resummation 
formula ((7j). 



Considering the quark-gluon channel, Eq. ([20]) can be recast in the following form: 



where we have used H.„n42„'Az) 



,DY{1) 
L qq^qq \ 



#f y(1) 5(1 - z) + 2 C^fc) (see Eq. ([11])). The relation flST]) can 
be used to determine Cqg (z) from the knowledge of Tif^gl and of the q F resummation coefficients 
at 0(as). Inserting the first-order results of Eqs. (jT4"]) - ([TE]) in Eq. ([3"T]) . we explicitly have 



Cf g \z) + \ H ^z{l-z) 



njDY(2) (A 
rL qq*^qg\ 4 ' ) 



zlog(z) + - (1 



+ 



^-4 
2 



z(l 



(32) 



where HfZlgl is given in Eq. (]2Tj) . Note that the right-hand side of Eq. ( |3T]) (or Eq. (132]) ' 



is resummat ion-scheme independent. Analogously to Eq. (TIB]) , the dependence of Cqg on the 
resummation scheme is thus parametrized by the first-order coefficient Hq Y on the left-hand 
side of Eq. ([32]). 



*The other non-vanishing entries are obtained by the symmetry relation 'Hr^l ab 
the second-order matrix W^^^l are vanishing because of Eq. 



H 



DY 



Several entries of 
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The process- independent coefficient functions Cf q \z), C yz L(z) C^Mz) and C^Uz) are obtained 



.(2), 



y(2), 



t(2), 



qq' 



<?<? 



analogously to Cqg(z). Considering the flavour diagonal quark-quark channel, Eq. ( T2"Uj) gives 
2Cf q \z) + 5{l-z) 



H° Y{2) -\{H? ni) ) 2 



+ - H DY{1) U DY(1) -(z) 

' 2 1 rL qq-^qq\ 4 'J 



rL qq<-qq\ /C ' ) ^ \ rL qq-^-qq w rL qq+-qq 



( z ) 



(33) 



where the right-hand side of Eq. ( 1551) is expressed in terms of resummation-scheme independent 
functions. Inserting Eqs. f|T4|) -flT6l) in Eq. (155j) . we explicitly obtain 



2Cj?(20+<5(l-s) 



; Cl - (-k 2 - 8) //" y w 



oy(2) r x _ o F 

qql-qqK^J j 



7T 2 -8) 2 



5(1 -z)- 



+ (tt 2 - 10) (1 - z) - (1 + z) Inz 



(34) 



where Tif^gg is given in Eq. (1251) . We observe that C q 2 q(z) includes a resummation-scheme 
dependent part that depends on Hq Y and H q Y{ ~ 2 \ We also recall [1] that the resummation- 
scheme invariance relates Cqq , H q Y ^ and the third-order coefficient of the quark form 
factor. 



Considering the flavour off-diagonal quark-quark channel in Eq. ( [20]) . we obtain 



Cqg( Z ) —Hqq*-qq( Z ) ) ^qq'{ Z ) ~ ^- qq Jql'( Z ) i @q q> i Z ) ~ ^qq^qq' ( z ) 



qq<~qq' 



.(2), 
'qq' 



qq^-qq' 



(35) 



where ft^X" 

.(2), 



ft 



DY(2) 

qq^qq > 



and are given in Eqs. (l24 l) -( [26l) . The off-diagonal second-order 



qq^qq 
.(2) 



coefficients C$j(z), C ( qq >{z) and C$ (z) are resummation-scheme independent. From Eq. (1261) we 

C q i{z). The equality between C qq i{z) and C q qi[z) is expected to 



observe that we have C^J,(z) 



be violated at higher perturbative orders (i.e., we expect C g ',(z) ^ C g l{z)). 



In this paper we have considered QCD radiative corrections to vector boson production in 
hadron-hadron collisions. We have presented the analytic result of the NNLO calculation of the 
vector boson cross section at small values of qx (see Eqs. d2J and (JSJ)). The NNLO result is 
compared (see Eq. (I22p ) with the predictions of transverse- momentum resummation. The com- 
parison gives a second-order crosscheck of the all-order resummation formula (jTJ), and it allows us 
to determine the previously unknown resummation coefficients at 0(a\). These are the coefficient 
functions H q Y ^ b {z) (see Eqs. CT-M) and the related coefficients Cf 9 \ C { qq \ C^,(z), C$j(z) 

and C y q >, \z) (see Eqs. (132|) . (1341) and (1351) ). which control the dependence on the rapidity of the 
vector boson. The knowledge of these second-order coefficients is relevant for phenomenological 
applications of both resummed and fixed-order QCD computations. These coefficients have been 
already implemented in resummed calculations of the inclusive qr distribution at full NNLL ac- 
curacy [22]. Using the method of Ref. [12] , the same coefficients have been used to perform the 
fully-exclusive NNLO perturbative calculations of Refs. [13] and 
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